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The momentum spectrum and the polarization of charginos and neutrahnos in squark decays are 
, affected by gluon radiation in the decay process q qx{g). Wc determine these corrections and 

^ , study their impact on the [q£] invariant mass distributions for leptonic x decays. The higher-order 

corrections, though small in general, can be sizeable near pronounced edges of the final-state distri- 
butions. 



Abstract 



1 Introduction 



The properties of supersymmetric particles, colored and non-colored, can be studied at the LHC in 
cascade decays of squarks and gluinos, which have large production cross sections in high-energy hadron 
collisions. In particular, squarks that are lighter than gluinos decay exclusively to quarks and charginos 
or neutralinos, q ^ qx- Subsequent decays of neutralinos into final states with two leptons are especially 
well suited to investigate properties of the supersymmetric spectrum. Invariant mass distributions of the 
final-state leptons and jets can be exploited, for example, for precision measurements of supersymmetric 
particle masses (see e.g. [1-4]). Moreover, correlations among the particles in the decay chains allow 
spin measurements of intermediate particles [5-8]. 

In the present article we concentrate on one aspect in the multitude of theoretical challenges involved 
in the description of supersymmetric cascade decays, i.e. the super-QCD (SQCD) radiative corrections 
at next-to-leading order (NLO) and their impact on energy distributions, polarization degrees, and on 
the shape of quark-lepton invariant mass distributions. Precise theoretical predictions for the shape of 
the distributions are required also at some distance away from kinematic edges and thresholds as to 
allow for sufficiently large, statistically significant event samples. Moreover, the shapes may be needed 
to resolve possible ambiguities in mass determinations from kinematic endpoints [9], they are crucial 
for the spin determination [5], and they allow to test models of supersymmetry breaking [10]. We 
investigate the impact of QCD radiation on the shape of distributions in NLO perturbation theory and 
derive, largely, analytical results which provide us with the proper understanding of the SQCD effects. 
Extending the Monte Carlo calculation of QCD effects in Ref. [11], we include gluon radiation off the 
final-state quarks. [Hadronic jet fragmentation which must finally be added is beyond the scope of the 
present study, but the effect can be estimated by means of QCD power corrections.] 

We shall, for definiteness, focus on cascades of L-type squarks, qi, of the first two generations 
without noticeable L/R mixing and consider the decay chain 



which has been discussed in great detail in the literature for the SPSla/a' benchmark points [12, 13]. 
The neutralino Xi is assumed to be the lightest supersymmetric particle (LSP) and thus is stable in 
R-parity conserving theories. Since the Born terms can be factored out, the NLO analysis applies, 
without modifications, to all neutralino and chargino decay chains of the same form. R-type squarks, 
(jR, decay in large parts of the supersymmetric parameter space, particularly at the benchmark points 
SPSla/a', to the LSP, qji — > qxi, and do not develop cascades. 

Spin correlations are important for the proper description of the cascade decay (1.1). Since only 
the gaugino components of the neutralino X2 couple effectively to the squark/quark current if quark 
masses are neglected, the qi squarks decay to left-handedly polarized q^ quarks. By angular momentum 
conservation also the neutralino X2 polarized left-handedly at the Born level. Gluon radiation in the 
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Figure 1: Generic Feynman diagrams for virtual and real super-QCD corrections to the decay q qx. 

decay process leads to an admixture of right-handed neutrahnos. The polarization is reflected in the 
energy and angular distributions of the decay products. 

The paper is organized as follows. In section [2] we shall present the results for the NLO SQCD 
corrections to the decay qi qx and their impact on the x energy distribution and polarization. The 
phenomenology of the complete qL decay chain (1.1), including the spin correlations, is discussed in 
section 13.11 In section 13.21 we analyze the decay chains in the broader framework of squark production 
at the LHC, including i2-squarks and anti-squarks. We conclude in section [H 

2 Super symmetric QCD corrections to the decay q — ^ qx 

Squark and gluino decays are affected by SQCD corrections. For their mutual decay modes, dominating 
if kinematically allowed, the NLO corrections were determined in Ref. [14]. The SQCD corrections to 
the partial widths for the electroweak squark decay channels. 



were calculated in Refs. [15-18]: Virtual gluon and gluino/squark exchanges renormalize the qqx vertex 
and additional gluons are radiated in the final state, cf. Fig. [TJ 

Vertex corrections neither change the x momentum nor its polarization state, but both are affected 
by gluon radiation. The x liiis spectrum of the two-body decay becomes a continuous spectrum (sec- 
tion 12. ip . Since in radiative decays the x is not back-to-back anymore with the polarized zero-mass 
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quark g, spin-flipped x states are admixed (section 12. 2p . While the effect on the polarization is small, 
the overall impact of gluon radiation on the final-state distributions, particularly near edges, is quite 
significant, nota bene in view of the envisaged high-precision analyses of these cascade modes at the 
LHC. 



2.1 X Energy distribution 



The Dalitz-plot distribution for the unpolarized decay q 
width Fb, 



qXdi normalized to the Born-level decay 
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(2.2) 



is expressed by Xq^g^^ = 2Eq^g^y^/Mq^ where Mq is the squark mass and Eg^g^^ are the energies of the quark, 
the gluon, and the x in the squark rest frame with Xy. + Xg + Xg = 2 and 2^/k < Xy < I + h. A useful 
abbreviation is the mass ratio k = M'^/M?. The partial width in Born approximation, Refs. [16,17], 



TB = -Mqil-Kff?^, 



(2.3) 



^--yv- ■ / Jqx 

depends on the chargino and neutralino mixing parameters [19,20] in the coefficients fqj^. Gluon radiation 
collinear with the final-state quark is described by the first term in Eq. (|2.2p since l + K — Xy = M^g/M?, 
where Mqg is the [qg] invariant mass. The second term damps the infrared singularity. 
Gluon radiation reduces the energy Xy. of the x according to the distribution 
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X,- 



4:K denotes the scaled neutralino momentum in the squark rest frame. 



Collinear and infrared gluon emission generate singularities at the maximum x energy Xy x^' 



1 -\- K. The leading contribution in 6 



Xy can be resummed by applying techniques analogous 



to thrust distributions [21-23] or lepto-quark decays [24]: 
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with the sub-leading coefficient 
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(2.6) 



Multiple gluon emission bends the distribution to zero at the kinematical boundary, corresponding to 
the Sudakov suppression of a final state without any gluons. Depending in detail on the effective value 
of Os, the turn-over point is very close to the maximum of Xy. 

The X energy distribution is displayed in Fig. for L-squark decays to x^^ with parameters 
adopted from SPSla' (i.e. k = 0.108) and the QCD coupling set to as(Mg) = 0.093. 
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2.2 X Polarization vector 

Gluon radiation at NLO also affects the x polarization. At the Born level the x polarization in decays 
must be left-handed. While the x spin does not flip for collinear and infrared gluon emission, because 
the intermediate left-chiral quark is effectively on-shell, the x spin can flip for non-collinear /non- infrared 
final-state configurations, with the angular momentum balanced by the emitted hard gluon. 

The polarization vector V[x] is defined in the x rest frame. The x-axis is chosen as the x flight 
direction in the q rest frame, the [x, ?/]-plane is identified with the decay plane, the y-axis pointing to 
the half-plane of shortest rotation from the x to the q momentum vector. The polarization vector lies 
within the decay plane, with components parallel, = Vx, and transverse, Vj_ = Vy, to the x flight 
direction while the perpendicular polarization component vanishes {Vz = 0). 

The two components of the polarization vector for fixed Xj^ can be cast into the form 



4 (2 - x^) {k [4 - x^] - x^) log (|-^^±^) + 



/AA, (2.7) 



^^(^x) = IjT^^T^^^— + '^^^ - 16Vl + «-^x(2 - ^x) - (32 - 3x^)] /M . 
o \L fi,; I -\- K — Xy^ 

The normalization factor J\f denotes the distribution {^^T ^)~^dV / dx-^ at the point x-^. For fixed 
Xy- < x'^^^ the polarization vector is independent of Og to leading order and depends only on the 
radiation dynamics. In the infrared and collinear limits of gluon radiation, — > x^^^, the Born values 
are approached again. 



-Pii -1 + -1 

" 2(l-K)2log(l/5) 

^ VI 

8 (1 - k) log(l/5) 



for x<^ ^ xf'''' = 1 + K . (2.8) 



In the opposite limit, x^ — > x™™ = 2-^/^, the polarization components approach the values "Pn and 
V± — > — 7r/4. The two components of the polarization vector are displayed in Fig. [2]J^b), again for the 
mass pattern of the SPSla' scenario {k = 0.108). 



2.3 Average x energy and polarization 

The values of the average energy and the polarization components of x for the inclusive NLO prediction 
are obtained by integrating over the Dalitz plot and including the 2-body decay limit. Both, the 
QCD corrected 2-body and the 3-body part of the calculation are infrared divergent. Analogously 
to the methods employed in Refs. [16, 17], small quark and gluon masses are introduced to regularize 
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Figure 2: (a) The differential decay rate for — > qgX2 function of the X2 energy Xy. at the 
benchmark point SPSla' (My./Mq = 0.329j. The curve labeled "resum" shows the damping near x™'*^ 
from multiple gluon radiation according to Eg. ( 12 .5)) . (b) The longitudinal and transverse polarization 
of the neutralino X2 ■ 



the singularities. It is sufficient to know the regularized form of the different contributions for the 
unpolarized case. Since the polarization effects are not infrared-sensitive, the only divergent integrals 
are those already encountered in the unpolarized results. The averaged observables (x^^), ("Pn), and {V±) 
globally reflect the deviation from the Born prediction to order a^, and they are given by 

and 

{Vu) = -1 - ^ — 7-^-^ \4.tt^{k - 3)k - 48{k - 3)kU2 (V^) + 12(1 - k[8 - 3k])U2{i^) 
" 9 vr (1 — Kj^ 

-3^/^ (22[1 + k] - Vk[45 - k]) +6(1 - k)(11 - 3k) log (\/^+ l) 

+3 (4[1 - k]2 log(l -k) + k[4 - 3k:]) log k] , (2.10) 

= ^7-^^[4 + 41og k-4k- (1 + V^) (15-k)E-2(k-7V^-8)K] . 

The function Li2 is the usual dilogarithm while K and E abbreviate the complete elliptic integrals of 
the first and second kinclll, respectively, evaluated for the argument [(1 — y/n) / (1 + y^)]^. 

For the reference point SPSla' (k = 0.108) adopted in the numerical analyses of the next section, 
^The integrals are defined as K{z) = /q"^^(1 — ^sin^ and E{z) = — 2; sin^ a;)^''^da::. 
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the mass-dependent coefficients of the 0{as) contributions take values 

(x^) = x^'^" - 0.32 — , (2.11) 

and 

= -1 + 0.02 — , {V±) = +0.08^. (2.12) 

TT O 

The NLO SQCD effect on the average longitudinal component of the x polarization vector, ('Pn), is very 
small. This follows from the fact that in the collinear and infrared regions, in which the density of the 
Dalitz plot is maximal, the deviation of the polarization vector from the Born approximation approaches 
zero. The average perpendicular polarization {V±) is positive, dominated by the contributions from the 
region of large x energies. 



3 Phenomenological results 
3.1 The Ql decay chain 

For parameter points like SPSla/la', the decays to lepton/slepton pairs of the first two generations, 

^o^i++i- and + (3.1) 

are ideal spin analyzers [R denoting right chirality] since the decays into the L-type state If^ are kine- 
matically forbidden, and L/R slepton mixing as well as the X2 higgsino component are suppressed. 
[Small corrections due to non-zero lepton masses can easily be included.] 

For the sake of clarity, we first focus on these decay modes of the neutralino. Since the positron 
is left-handedly polarized, it will be emitted preferentially anti-parallel to the x spin direction, 

1 r 1 

t;-^ ^ = -(l-P[x]cos0,), (3.2) 

i a cosPg 2 

where 9s is the angle between the x spin vector and the 3- momentum, and Vlx] the x degree of 

polarization. The angular correlation Eq. (j3.2p gives rise to an increase of the invariant mass of the 
[qi^] pair in the decay chain 

qL^q + xl + {9)^q + e^ + + {9) (3.3) 

compared with the isotropic distribution. [In this cascade the lepton I is generally termed "near lepton" 
as opposed to the "far lepton" emitted in the subsequent slepton decay I ^Xi; which will be included 
later.] 

We have implemented the NLO results presented above in a flexible parton level Monte-Carlo pro- 
gram to calculate arbitrary differential distributions. Both phase-space slicing and subtraction methods 
have been employed to allow for internal cross checks. 
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In the phase-space sHcing approach, approximate real radiation matrix elements are integrated an- 
alytically in the soft/collinear regions using mass regularization. After adding the mass-regularized 
virtual corrections, the result is finite and the regulating quark and gluon masses can be set to zero. 
The hard gluon region is integrated by standard Monte-Carlo techniques. 

Alternatively, a subtraction method has been applied which is particularly suited to construct NLO 
parton-level Monte-Carlo programs. In the present simple case the entire 3-particle differential decay 
width can be used as subtraction term. For the observables O which we consider, the expression in 
NLO can be cast in the form 

(O)nlo = / dtiPs -02]+ [ [dt^ + [ dr§]02 . (3.4) 

Here, dV^ is the differential 3-parton decay width, dT2 the vertex-corrected 2-parton decay width, both 
normalized to the NLO width F. Moreover, On denotes the observable O calculated for either n = 3 or 
n = 2 parton final states. The projection from the n+1 onto the n particle phase space is straightforward: 
at any given 3-body phase-space point, O2 is evaluated for a neutralino momentum which shares the 
direction of flight with the neutralino momentum of the 3-body phase-space point. The simple 2- 
body kinematics fixes the neutralino and quark momenta completely so that O2 can be calculated 
unambiguously. The first integral in Eq. (j3.4p is finite by construction and can be integrated by Monte- 
Carlo techniques. The second and the third integral can be evaluated analytically [in mass-regulated 
form to isolate the infrared/collinear divergences, which cancel in the sum of the virtual contribution 
and the integrated subtraction term]. Since polarization effects are not infrared sensitive, it is sufficient 
to use the unpolarized differential decay width dT^ as a subtraction term even for observables including 
spin correlations. In this case, O2 is evaluated with the tree-level polarization vector, which provides 
the correct limit of the 3-body kinematics in the soft/collinear regions. The results derived by using the 
subtraction method agree with those obtained from phase-space slicing. 

When gluon radiation is switched on, the [qi^] invariant mass distribution is softened. To quantify 
this effect, the (/-jet must be defined properly in an infrared-safe way. We will combine the q and g partons 
in the collinear and infrared regimes to a single ^-jet if their scaled invariant mass y = M'^[q = qg]/M? 
is less than yc = 0.01. For events with non-collinear hard gluon emission we identify q with the leading 
parton jet, i.e. the jet with the highest energy in the squark rest frame, either q or, rarely, g. 

The invariant mass distributions M'^[q£^] are displayed as solid curves in Fig. [3j In the upper left 
panel the distribution of the positively charged "near" leptons, produced in ^2 decays, is shown. The 
dashed curve represents the invariant mass distribution at the Born level for comparison. The peak 
of the distribution at high invariant mass is more pronounced than the triangular-shaped tree-level 
distribution with spin correlations ignored [25]. As expected, the NLO corrections lead to the transfer 
of events to smaller invariant masses compared to the tree-level prediction. 

The transition from the 1-jet to the 2-jet class of events gives rise to the tiny kink in the [^-^near] 
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Figure 3: The [ql^] invariant mass distributions for positively (left panels) and negatively (right 
panels) charged leptons is shown at SPSla' [specifically for ul masses and first and second generation 
leptons with vanishing masses]. The upper panels include only "near" leptons while the lower panels 
display the sum of "near" and "far" leptons. The NLO predictions are represented by the solid lines 
for y-cut parameters yc = 0.01 and 0.002, the Born result [approximately equivalent to NLO with 
2/c ^ 0.05/ hy the dashed line. 
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distribution in Fig. [3l Varying the experimentally defined jet-resolution parameter yc reproduces anti- 
cipated modifications of the mass distributions: the distribution is softened for smaller yc, and the peak 
rounded off, exemplified in Fig. [3]for yc = 0.002, as a larger fraction of events with respect to the default 
choice yc = 0.01 moves from the 1-jet to the 2-jet class. For larger yc the distribution is sharpened, 
moving back towards the Born distribution, as more events are attributed to the 1-jet class. The NLO 
distribution becomes almost identical to the LO prediction for yc > 0.05. Quite similar results for the 
invariant mass distributions are found when using a cone algorithm [26] defined in the squark rest frame. 

Parallel to the x energy distribution, exponential Sudakov damping of the invariant [q£^] mass 
distribution is expected in a margin of one percent below the pronounced kinematical edge when multiple 
gluon radiation is switched on. Likewise, hadronic jet fragmentation will give rise to smearing effects, 
generally accounted for, see e.g. Ref. [27], by shifting the perturbative prediction of the distribution 
downwards by a relative amount of approximately 1 GeV /Mq. For squarks decaying through electroweak 
channels, as at the present reference point, non-zero width effects oc Tq/Mq are also limited to about 
one percent, on both sides of .he kinematical edge, though! These effects will round off the peak of the 
distribution at the edge for yc > 0.01, with the overall impact limited however to a margin at the one- 
percent level for squark masses in excess of 500 GeV. With details depending on the jet definition, the 
modification will qualitatively be similar to the thrust distribution in jet analyses of e^e~ annihilation 
at high energies [27]. The explicit Monte Carlo simulation of QCD effects in Ref. [11] includes only 
gluon radiation off the squarks which turns out, as expected, to have a very small effect on the shape 
of the distributions. 

The experimentally observed lepton l'^ may also be the decay product of the slepton in the chairj^ 
X2 ~^ ^ij^R ~^ ^R^R Xi- Foi' this "far" lepton l'^ , the daughter of a spinless parent, the spin character- 
istics in the distributions are largely washed out and only the average energy reduction of the slepton 
due to spin correlations in the primary neutralino decay is effective. 

Adding up the contributions from "near" and "far" leptons, which cannot be separated exper- 
imentally on an event-by-event basis, we arrive at the distributions displayed in the lower left panel 
of Fig. [3l Due to the rather small mass difference between the selectron/smuon and the LSP Xi i^i 
SPSla', the lepton from the slepton decay is relatively low-energetic. Thus the "far" lepton adds to the 
distribution only at low invariant masses and the signature of the "near" lepton is still significant at 
high invariant masses. 

If, instead, negatively charged leptons are observed the shapes of the two corresponding distributions 
are rather different. Since the negatively charged leptons are right-handedly polarized the sign in 
Eq. ()3.2p flips and the "near" £~ is preferentially emitted parallel to the ^-jet, in contrast to the anti- 
parallel emission of positively charged leptons. The energy of the negatively charged lepton l~ is 

^The accuracy of the narrow- width approximation for cascade decays has been addressed in Ref. [28]. 
^Electroweak radiative corrections to the decay X2 l^l^Xi have been presented in Ref. [29]. 
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Figure 4: The asymmetry Eq. 113. 5\) as a function of the [ql^] invariant mass at SPSla! for q^ 
decays [specifically for ul masses] into first and second generation lepton pairs (no anti-squark decays 
included). 



therefore, on average, reduced in the squark rest frame for polarized X2 decays. For both reasons, the 
[qi~] invariant mass is reduced on average compared to the prediction without spin correlations [25]. 
The reduction is slightly enhanced by gluon emission as demonstrated in the upper right panel of Fig. [3l 
The distribution for the sum of "near" and "far" i~ leptons is shown in the lower right panel. Its shape 
can easily be distinguished from the results at large invariant masses. 

Spin correlations are clearly important for the description of the invariant mass distributions; without 
taking spin correlations into account, all panels in Fig. [3] would be identical. Comparing the final 
distributions in the lower panels an asymmetry between positively and negatively charged leptons, 

is predicted indeed [5]. [The asymmetry is washed out somewhat in pp collisions at the LHC where 
also anti-squarks q* are produced; see section [3.2] , Fig. U] shows the parton-level asymmetry for the 
qL decay chain in SPSla'. In particular for high invariant masses the dominance of positively charged 
leptons is evident. The NLO SQCD radiative corrections included in our decay analysis have a very 
small influence on the asymmetry. This is evident from the lower panel of Fig. [Sj near the kinematic 
endpoint where the corrections to the invariant mass distribution of the positively charged lepton are 
biggest, the asymmetry is close to unity and changes are marginal as a result. 
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3.2 Squark cascade decays in pp collisions at the LHC 

The more detailed numerical analysis presented in this final section is also based on the reference 
point SPSla' [12] in the mSUGRA scenario. This point is located in the transition zone from the 
bulk to the coannihilation region in the [Mx/2)-^o] gaugino/scalar mass plane. It corresponds to the 
region of maximum probability in the analysis of electroweak precision observables [30,31], driven by 
the contribution to the muon anomalous magnetic moment, (7^ — 2; the reference point is compatible 
with the WMAP measurement of the density of cold dark matter. The parameters of this mSUGRA 
scenario are defined as follows: gaugino mass Mx/2 = 250 GeV, scalar mass Mq = 70 GeV, and trilinear 
coupling Aq = —300 GeV at the unification scale; moreover tan/3 = 10 and sign(/x) = +. After 
evolution down to the Terascale the squark, gluino, slepton, and neutralino masses that will be used in 
the phenomenological analysis are collected in Tabled) While Suspect [32] has been employed for the 
computation of the entries in the table, the Spheno [33] and SDFTSUSY [34] spectra are quite similar, 
with differences at the level of 1% and less. 

Table [1] also presents the branching ratios which are based on NLO calculations as implemented 
in Susyhit [35] / Sdecay [36]. The L-squarks decay predominantly into the charginos Xi ^^'^ the 
neutralino X2- -R-squarks decay almost exclusively into the x? LSP which is bino-like. Hence, almost 
all of the squarks decaying to X2 ^^^^ indeed L-type squarks. Only a small admixture of i?-type squarks 
diminishes the observed spin correlation and the asymmetry, Eq. (|3.5|) . generated by the X2 decay chain. 
More than half (i.e. 57%) of the neutralinos X2 decay visibly into charged slepton-lepton-pairs. Since 
left-handed sleptons of the first two generations are heavy, the decays are dominated (BR = 0.527) by 
decays into mixed fi states [37] with a reduced mass eigen-value, but still a significant fraction decays 
to e/j,/iR (BR = 0.024 for each mode). If the e,/i leptonic decays are used as experimental signals, the 
transverse momenta of the cascading r — > e, /i final states are strongly reduced. A large fraction of the 
T ^ e, decay events will thus be lost in experimental cuts, about 90% for a transverse momentum cut 
of 10 GeV. 

Since experiments will not distinguish quark from anti-quark jets, experimentally observable distribu- 
tions will also receive contributions from anti-squark decay chains. For anti-squarks, by CP-invar iance, 
the role of positively and negatively charged leptons is interchanged and the asymmetry A is reduced as 
a result. At the Tevatron, the equal rate for squark and anti-squark production even drives A to zero. 
However, the (non-CP-invariant) proton-proton initial state at the LHC ensures that more squarks than 
anti-squarks are produced and a non-zero asymmetry is predicted. 

At the LHC, squarks are produced either directly in quark/anti-quark/gluon collisions, or they are 
decay products of gluinos as Mg > Adq for the SPSla' reference point. The production cross sections 
for both species have been determined including the next-to-leading order corrections of Ref. [38] as 
implemented in Prospino [39]. The NLO calculation of Ref. [38] has been performed for mass-degenerate 
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Table 1: Masses and decay branching ratios most relevant for evaluating the SPSla' squark/gluino 
cascades. 

L and i?-squarks and mass-degenerate first and second generation squarks of up and down type. Since 
the mass differences are in general small, in particular for the spectrum of the SPSla' reference point 
where [M£ — M^]/[M^-|-M^] ~ 4%, this is a valid approximation. In our tree-level calculation, to which 
we apply the K-factors, the exact kinematics is used. For the parton densities the CTEQ5 set [40] has 
been adopted and the renormalization and factorization scales have been set to the average mass of the 
produced particles. The production cross sections are displayed for squarks and anti-squarks of the first 
two generations in Table [2j The pp initial state at the LHC clearly favours the direct production of 
squarks over the direct production of anti-squarks. Though equal numbers of squarks and anti-squarks 
are generated by gluino decays, the ratio of the number of anti-squarks over squarks is only slightly 
diminished because the branching ratio of gluinos into first and second generation L-squarks is small. 

Table[2]also includes the cross sections times branching ratios a{q) x BR(g — > X2) ^id) ^ BR(^ — > 
Q ^ X2) neutralino X2 production from squark decays, and the final cross sections times branching 
ratios for the signal cr[e,/i], i.e. lepton production from squark cascades, with r decays separated. With 
a nominal integrated luminosity of 300 fb~^, a sample of roughly 2 • 10^ events is accumulated for the 
SPSla' benchmark for e, /x in the first two generations and, before experimental cuts, an initial number 
of 1.5 • 10^ r events decaying leptonically into either electrons or muons. Although these calculations 
are of purely theoretical nature, they provide nevertheless a solid platform for estimating expectations 
before experimental simulations will finally include selection cuts, efficiencies etc. 

Using these estimates of production cross sections and branching ratios, the invariant mass distri- 
butions for the SPSla' decay chain can be calculated. Fig. [5] displays the LO and NLO [jet^^] invariant 
mass spectrum for positively and negatively charged leptons as predicted for the LHC. The distribu- 
tions are normalized to the signal cross sections and are shown separately for X2 decays to first /second 
generation leptons / sleptons and for r /f decaying leptonically. L/R mixing of the r sleptons and the X2 
higgsino component exacerbate the analysis of the r/f channel considerably; all these effects have been 
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Table 2: Cross sections for the production of squarks of the first two generations and gluino production 
in pp collisions at the LHC, and the production cross sections for the neutralino X2 cascades. 
Also shown are the cross sections for R-sleptons after summing up all cascade channels. The resulting 
cross sections for electron/muon pairs from direct electron/muon production in the decay chain and 
from tau decays are collected in the right-most column. 

taken into account in the numerical analysis. However, the leptonic r decay signal is strongly reduced to 
a level of 10% if a transverse momentum cut of 10 GeV is applied as indicated by the dotted lines. [For 
simplicity, the cut is applied in the squark rest frame. This approximation is adequate since squarks are 
produced predominantly at small velocities.] The experimentally observable distributions are clearly 
distinct despite of the anti-squark admixture to the squark sample. As a result of the anti-squark ad- 
mixture, the lt~ asymmetry (j3.5p is slightly reduced at the LHC compared with Figured for pure 
qi, decays, with marginal NLO corrections as before. However, the asymmetry is still quite significant 
for large [jet£] invariant masses as demonstrated in Refs. [5-7], even after including detector effects. 

4 Conclusions 

Edges/thresholds, distributions and correlations of particles in cascades play an important role in de- 
termining properties of supersymmetric particles, i.e. masses and spins. Lepton-lepton and jet-lepton 
invariant mass distributions are particularly useful instruments in this context. 

In the present article we have studied how the parton final states in g — > decays are affected by 
super-QCD corrections at next-to-leading order, in particular by the radiation of gluons in the decay 
process. The corrections modify the energy as well as the polarization of the particles. By performing 
an analytic next-to-leading order analysis we have provided a detailed understanding of origin and 
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Figure 5: Invariant [ql] mass distributions, yc = 0.01, for muons and electrons produced directly in 
the neutralino and slepton decay. For muon and electron pairs produced from tau decays in the decay 
chain, the reduction of the NLO cross sections for a transverse momentum cut of 10 GeV is indicated 
by the dotted lines. 

size of these effects. Moreover, we have studied particular decay channels and [jet^^] invariant mass 
distributions, exemplified for the SPSla' benchmark scenario, in — > qi~^i~ Xi decays. We have found 
significant gluonic corrections to the shape of the invariant [jet i^] mass distribution near the pronounced 
edge at the kinematic endpoint. 
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